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Abstract 

We present a complete order a s analysis of the process pp — > 77 + X with 
polarized initial states, previously studied by us with leading order structure 
functions. We include in our calculation new sets of parton densities evolved in 
NLO QCD, such as those of Gehrmann and Stirling and of Gluck et al., which 
incorporate the new anomalous dimensions of Mertig and Van Neerven in the 
evolution equation. A detailed phenomenological analysis is also given which 
includes the photon isolation. Our results indicate that the asymmetries, 
although not very large, should be substantial enough at RHIC energies to be 
measurable in future planned experiments. 
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I. INTRODUCTION 



The study of the spin structure of the nucleons is a fascinating subject which has at- 
tracted a lot of attention both from the theoretical and from the experimental side in recent 
years. From a phenomenological viewpoint, for instance, the EMC result has called for a re- 
analysis of the connections between the traditional parton model and the QCD description 
of polarized DIS, based on leading-power factorization theorems and on applications of the 
Operator Product Expansion. 

There seems to be more agreement now in the literature on the interpretation of the 
EMC and of other more recent results than before, and it seems obvious that even in a 
perturbative framework (such as the OPE), subtle issues connected to the renormalization 
of anomalous operators, their scale dependence etc., had been overlooked in the past. Sum 
rules, largely inspired by the naive quark model, have also provided a bridge between the 
theory and the experiments. 

It is expected in few years time that the information gathered from DIS studies of the 
polarized parton distributions will be supplemented and even extended by new results ob- 
tained from proton-proton colliders, such as the BNL Relativistic Heavy Ion Collider, RHIC. 
In particular, in these experiments a direct gluon coupling will allow direct measurements 
of the contributions of the gluons to the spin of the nucleon. 

Recently, new sets of polarized parton distributions have been generated by various 
groups which incorporate the effects of the evolution equation up to 0(a s ). This has been 
possible thanks to the ground-breaking calculation of Mertig and Van Neerven, later ver- 
ified by Vogelsang || who have derived the 0(a s ) radiative corrections to the anomalous 
dimensions of the evolution. 

We have presented a study of the process pj)->77 + Xto next-to-leading order 



containing a brief phenomenological analysis of our results [17]. In this work we are going 
to extend this analysis for the same process by incorporating, for the first time, the new 
NLO sets of structure functions. Here most of our discussion, therefore, will be purely 
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phenomenological and we refer to our previous study for all the technical details and for a 
complete presentation of the methods involved in it. 

This work is organized as follows. In section II, in order to make our discussion self 
contained, we briefly review the structure of the evolution equations for the structure func- 
tions to order a s , and discuss in general terms the various parametrizations. We then move 
to a detailed phenomenological study of double prompt photon production, and examine 
in particular whether it will be possible to get any new information on the polarized gluon 
distributions AG from a study of this process at RHIC. We will also discuss the uncertainties 
and the scale dependence of NLO results, and we make predictions of the large transverse 
momentum (jpr) behaviour of the cross section at RHIC. Our conclusions are presented in 
section IV. 



II. THE NLO EVOLUTION 

Different sets of polarized structure functions to NLO have been proposed recently by 
Gehrmann and Sterling (G-S) |7J and by Gliick et al. (GRSV) 0. For the G-S sets, three 
different parametrizations are considered, while for GRSV two possible scenarios have been 
analyzed. In order to render our treatment self contained, we briefly summarize these new 
results. We start from the usual framework of a DIS (longitudinal) polarized scattering. 

In the QCD parton model we define 

q = g T + q l 

Aq = q^-q l (2.1) 

to be the unpolarized and polarized distributions respectively, q^ and q^ are the distributions 
which describe the probability for finding quarks with spin parallel or antiparallel to the 
longitudinally polarized nucleon. For instance, if we use values for the factorization scale 
and the renormalization scale both equal to Q 2 , we can write down the relation between 
the polarized structure function g\ and the polarized quark (Ag(x)) and gluon distributions 
(AG(x)) by 
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9l (x,Q 2 ) = f 1 ^ (C^x/y,a s (t))AE(y,t) + C? s (x/y,a s (t))Aq NS (y,t) 
Jo V 



>o y 

+ 2n f C g (x/y,a s (t))AG(y,t)), 



(2.2) 



where 



n f 

AE(x, t) = £ ( A ^(*, t) + Aq t (x, t)) 
i=i 

A <?^ = E ^'tt ( A <^, *) + 0) • (2-3) 

i=i \ e / 

are respectively the singlet and the non-singlet quark distributions, n/ is the number of 
flavors (i), each of charge e*, (e 2 ) = I]e 2 /n/ an d * — l°g(Q 2 /^- 2 ) 
The singlet part of the Altarelli-Parisi evolution is given by 

d AE(x, t) = /' ^ [P*(a;/s/, a s (t))AT,(y, t) + 2n f P qg (x/y, a s (t)AG(y, t) 



dt 2ix Jx y 

±AG(x,t) = *®[ 1 *l 
dt v ' 2tc Jx y 

while the non-singlet part evolves independently 



d AG(x, t) = ^M f'^l \ P S( x /y : a s (t))AY,(y, t) + 2n f P gg (x/y, a s (t)AG(y, t)\ . (2.4) 

Z7T Jx V L J 



At NLO the expansion in power of a s of the coefficient functions (C) and of the splitting 
functions (P) reads 

C(x, a s ) = C< >(x) + %-CP\x) + 0(a 2 s ), 

In 



P(x,a s ) = P(°\x) + — pW(x) + 0(a 2 ). 

2% 



(2.6) 



The first moment of gi(x, Q 2 ) measures the expectation value of a combination of octet 
and singlet axial vector currents 

1 111 

gi(x)dx = — a 3 + — a 8 + -a , (2.7) 
o 12 do 9 



where 



a 8 = 3F - D 

= Ah + An + Ad + Ad - 2(As + As) = 0.579 ± 0.025; 



(2.8) 



a 3 = g A = F + D 

= Au + Au - Ad - Ad = 1.2573 ± 0.0028, 



(2.9) 



a = ^ A g + Ag = Ag 8 + 3(As + As). 

<i 

(2.10) 

Aq denotes the first moment of the corresponding distributions Aq(x). Under suitable 
assumptions ||, an analysis done in 1983 gives for the matrix elements controlling the (3- 
decay of the spin 1/2 hyperons 

F = 0.477 ± 0.012 D = 0.756 ± 0.01. (2.11) 

In a more recent analysis (1988) [|J 

F = 0.46 ±0.01 D = 0.79 ± 0.001, (2.12) 

which implies a = 0.06 ± 0.12 ± 0.17 if the EMC data (r£[< Q 2 >= 10.7] = 0.128 ± 0.013 ± 
0.019) or the SMC data (r?[< Q 2 >= 10] = 0.136 ± 0.011 ± 0.011) for the first moment 
of gi are used. This value is smaller than the value ao = 0.188 ± 0.004 expected on the 
basis of the Ellis- Jaffe sum rule. The values of a 3 and as are constraints equations which 
have been used in most of the LO analysis of the polarized structure functions done so far. 
These two constraints are obtained under different assumptions. For instance (|2.9| ) requires 
SU(2)f symmetry between the matrix elements of the charged and neutral axial currents, 
while (|2.8p requires an SU(3)f symmetry between matrix elements of hyperon decays of 
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charged and neutral weak axial currents. This last assumption has been criticized by Lipkin 
fl~3| . He has argued that the /3 decay of hyperons only fix the total helicity of the valence 
quarks. Therefore eqs. ( 2.8|) and ( |2.9| ), in this scenario, should be replaced by 



a 8 = Au v + Ad v 



(2.13) 



and 



0,3 = Au v — Ad v 



(2.14) 



respectively. These two scenarios, the standard scenario and the valence one, generate two 
different sets of parton distributions GRSV1 and GRSV2 in ref. ||. The expression of these 
two parametrizations can be found in the original work. Notice that these two scenarios, of 
course, require different assumptions on the contribution to the polarization coming from 
the sea quarks. 

In fact, the two expressions for the first moment of gi, assuming the validity of the 
Bjorken sum rule for (eq. |2.9| ) and of SU(3)f flavour symmetry are given by, 



1(F + D) + ^(3F - D) + i(Aa + As) 



« S (Q 2 



7T 



(2.15) 



for the standard scenario, and 

r? = 



1(F + D) + ^(3F -D) + ^(lOAg + As + As) 



« S (Q 2 



7T 



(2.16) 



for the valence scenario. Notice that in the standard scenario (eq. |2.15|) , a negative As has 
to be required in order to bring down to the experimental EMC or SMC result. As we 
just mentioned, the condition As = As = is in fact still too high to match the data 

In the valence scenario of Lipkin (eq. |2.16| ) - even with zero strange quark contribution 
- a negative Aq is needed in order to obtain the same reduction. In general, a finite As(Q 2 ) 
is generated by the evolution equations due to the non vanishing of the NLO anomalous 
dimension A'jgg. This happens for Q 2 > /1 2 , where fi is the renormalization scale to NLO 
where the evolution starts. In GRSV /i 2 = 0.34 GeV 2 . 



Let's remark, at this point, that our understanding of higher twist effects and renormalon 
effects at small Q 2 , which influence the sum rules and might affect the evolution, are not yet 
under control from a theoretical perspective. It is estimated, though, that at Q 2 = 2 Gev 2 , 



for instance, higher twist contributions should be of the order of f 0% |]T4 . 

It should be pointed out that it is quite common in the literature to find different ex- 
pressions for depending over whether the anomalous contribution is included or not in 
the definition of ao- Notice that this amount to a redefinition of AS of the form 

AS ^ AS - ^^-AG (2.17) 
2n 

which modifies the contribution to g\ which can be attributed to the net helicity of the 
quarks by an anomalous gluonic part. However, from a practical viewpoint, it is convenient 
to remove it from the LO expression of Ti and let it reappear through the evolution equation 



(in the MS scheme), now known to order a s , as an additional term in the splitting functions. 



The calculation of the NLO anomalous dimensions of ref. |J are performed in the MS 
scheme. 

Generically, the ansatz chosen by GRSV for the polarized parton distributions is of the 
form 

Aq(x,Q 2 ) ~x°(l -xfq(x,Q 2 ) (2.18) 
and a similar form for the gluon content, where the unpolarized distributions q(x, Q 2 ) are 



taken from [11|. The gluon distribution is only weakly constrained since, in LO, it doesn't 



appear in the expression of the asymmetries. Various arguments, based on Regge theory or, 



as in ref. [[L5[], arguments based on coherence effects, allow to estimate a behaviour of the 
form 

AC 

(2.19) 



G 

Notice that the extrapolation in x of the various distribution is a remarkable source of 
complexity in spin physics, given the fact that the measurements are performed at different 
values of Q 2 at each bin x and an intermediate AP evolution to common Q 2 is required. 



The G-S distributions in ref. are generated by global fits of the form 

xAQ = r]Ax a {l-x) b {l + -fx + pv^) (2.20) 

where three different parametrizations are given (here denoted G-S a, G-S b, G-S c). They 
differ in the parametrization adopted for the polarized gluon distributions. Specifically 

G-S a 7g = p G = 0, 
G-S & 7G = -1 Pg = 2, 

G-S c 7g = p G = -3. (2.21) 

We refer the reader to the original work for further details. In the next section we 
are going to discuss the implications of these parametrizations for the measurements of the 
double photon cross section at RHIC. 



III. NUMERICAL RESULTS 

We now turn to the numerical results for polarized isolated double photon production 
at RHIC. It is still not clear what isolation parameters will be used in the experiments, 
but preliminary studies |§ seem to indicate that an energy resolution parameter of e = 
0.5 GeV/py will be possible. We assume a cone size of R = 0.5, where R is defined as the 
radius of a cone in the pseudorapidity-azimuthal angle plane via the relation 



R = y (Ay) 2 + (A0) 2 . (3.1) 

If a parton with energy fraction greater than e times that of the photon falls into the cone, 
the event is rejected. 

As discussed in section II, we make use of the recent polarized proton distributions 
evolved in NLO by Gehrmann and Stirling H and Gliick et al [H. For the unpolarized proton 



distributions, we use the CTEQ3M set pi after checking that the GRV |ITJ distributions 



give very similar results. We also use fragmentation functions for the photon evolved in NLO 
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from ref. ||12|| . The value of Aqcd used is chosen to correspond with the parton distribution 
set used. For the electromagnetic coupling constant we use a em = 1/137 and we use the 
two- loop expression for a s . Unless otherwise stated we set all renormalization/factorization 
scales to /z 2 = ((pti) 2 + {Pt-i) 2 ) / 2 > where p^t and the transverse momenta of the first 

and second photons respectively. 

In Figs. la and lb we compare the different parametrizations of AG(x, Q 2 ) as a function of 
x at Q 2 = 100 GeV 2 . The two GRSV gluons are practically identical so we would not expect 
very different predictions for any cross section sensitive to AG from these two parametriza- 
tions. The three G-S distributions on the other hand, show much larger differences, and are 
thus likely to give different predictions for processes depending on AG. 

At the moment, it is proposed that RHIC will run at various centre-of-mass energies 
between v^S* = 50 to 500 GeV, thus we perform calculations for two values, = 200 
and 500 GeV. It turns out that the NLO parametrizations for the polarized distributions 
do not include charm quark distributions, thus in our study we neglect the contribution for 
charm quarks in both the polarized and unpolarized cases. It is expected that at RHIC cms 
energies, the charm contribution will not be very large. 

In keeping with the convention started by the WA-70 collaboration [ lEfl , and which has 



been used in all subsequent analyses, we place asymmetric cuts on the two photon p' T s. We 
require that the photon with the highest transverse momentum has pt > 10 GeV while 
the other is required to have pt > 9 GeV. This asymmetric cut ensures that the two-body 
contributions to the cross section, such as the Born contribution, are not favored over the 
three-body contributions. This differs from the way we calculated the cross section in ref. 



17fl , where we looked at the pr distribution of one of the photons while placing cuts on the 
other. In the present way of calculating the cross section, each photon which passes the 
cut contributes to the cross section. This means we have two entries for each event. Latest 
indications are that the maximum rapidity coverage achievable at RHIC will be in the range 
— 2 < y < 2, we thus restrict ourselves to this range. It is also suggested that they will not 
be able to reliably detect photons with pt < 10 GeV. 
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For comparison, in Fig. 2a we show the px distribution for the unpolarized case as pre- 
dicted using the GRV and CTEQ3M parton distributions at ^fS = 200 GeV. They clearly 
give very similar results in the region tested, hence, in all subsequent discussions we will use 
only the CTEQ3M distributions. In Figs. 2b and 2c we display the corresponding polarized 
cross sections. The results indicate that it will be difficult to measure this cross section 
beyond about pr = 20 GeV, even with the planned high luminosities at RHIC. More inter- 
estingly, Fig.2b suggests, as we expected, that GRSV distributions give very similar results 
and cannot be separated by this process. Fig2c suggests, on the other hand, that the three 
G-S distributions give substantially different predictions, and will probably be distinguish- 
able. The G-S a parametrization gives similar results to the two GRSV versions, as may be 
expected from Figla. 

In Figs. 3a and 3b we show the longitudinal asymmetries as predicted by the various 
parametrizations. The asymmetry, All, is defined by the relation 

All = j^, (3.2) 

the ratio of the polarized to the unpolarized cross section, and gives a measure of the 
spin dependence or spin sensitivity of the process. Again the two GRSV distributions, 
as expected, give similar results, while the G-S ones give clearly distinguishable results. 
The predictions also indicate that the asymmetries are not very large in the measurable 
region of the cross section, varying between 5 and 8% for the GRSV and 4 and 6% for 
the G-S a parametrization. It will require very high statistics measurements to measure 
this asymmetry, but this may be achievable at RHIC, as long as the cross section can be 
measured, given the planned detectors. 

At = 500 GeV the situation improves somewhat with regards to the size of the cross 
section. Fig. 4a compares the p T distributions at v^S* = 500 and 200 GeV for the unpolarized 
cross section. As expected, at higher cms energies the cross section is substantially larger. In 
Fig. 4b we show the polarized cross section at \/~S = 500 GeV. On the figure we also show the 
full prediction given by the G-S a distribution as well as the contribution form qg initiated 

10 



subprocesses. In the case of G-S a the contribution for the qq initiated subprocess (not 
shown) is negative. Corresponding distributions are shown for the GRSV 1 parametrization. 
In this case the qq process gives a positive contribution but it is substantially smaller than 
the qg one. The G-S a and GRSV 1 distributions predict similar cross sections, but the 
relative importance of the subprocesses is clearly very different. The most interesting aspect 
of this result from the point of view of sensitivity to AG is the fact that in both cases, the qg 
initiated process dominates. We should mention that in this calculation we do not include 
contributions for the higher order process (0(a%)) gg — > 77, preferring to keep consistently 
to 0(a s ). 

Figs. 4c and 4d show the asymmetries as predicted by the G-S a and GRSV polarized 
distributions. Included also are the contributions to the asymmetries from the qg and qq 
initiated processes. In Fig.4d we also include the asymmetry in LO where we have used the 
LO counterpart of the GRSV 1 parametrization. The LO cross section consists only of the 
process qq — > 77 plus the fragmentation processes qq — > gj plus qg — > qj. It can be said that 
the asymmetry is fairly stable under the higher order corrections, but the corrections are 
clearly significant. As we expected from the results shown in Fig.4b, the NLO asymmetry 
is dominated by the qg initiated process, this is the case when either distributions is used. 
On the other hand, the size of the asymmetry is not very different from that at = 200 
GeV, varying between 3 and 8% over the measurable range. 

As an indication of the stability of our predictions we show the scale dependence of 
the unpolarized cross section in Fig. 5. The cuts and distributions are the same as for the 
solid curve in Fig. 4a. By varying all renormalization and factorization scales in the range 
1/2 < n < 2, where \x = ^((Pti) 2 + G°T2) 2 )/2, we can vary the cross section by as much as 
20% in the region 10 < p? < 35. This is still a substantial uncertainty in the predictions, 
but it is a significant reduction over the corresponding figure of 60% for the LO predictions. 

Finally in Fig. 6 we show the K-factor as a function of pr where, 

dcr(LO) 



da (NLO) 
dpjp 
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at v^S* = 500 GeV, for the unpolarized cross section. This confirms that the NLO corrections 
are indeed quite significant at the pr values relevant at RHIC. 

IV. CONCLUSIONS 

An 0(a s ) NLO calculation of the process pp — > 77 + X was presented, where for the 
first time polarized parton distributions evolved in NLO QCD is used. The two photons 
were also isolated using plausible isolation parameters. We found that the cross section for 
the precess will not be very large at RHIC, particularly at the lower CMS energies, but 
it will nevertheless still be measurable in the lower p\ region. The new polarized parton 
distributions predict an asymmetry of between 3 and 8% in the measurable region for the 
process. While the asymmetry cannot be said to be large, given sufficiently good statistics, 
it should still me measurable. Given this possibility, we found that a discrimination between 
extreme parametrizations of the polarized gluon distribution AG should be possible. This 
is due in particular to the dominance of the subprocess qg — > 775 over the qq annihilation 
process. We thus conclude that this process, if studied at RHIC could prove useful for 
supplementing information on the polarized distributions obtained form other sources. 
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Figure Captions 



[1] (a) The polarized gluon distribution in NLO as a function of x at Q 2 = 100 GeV 2 
for the G-S a (solid line) GRSV1 (dashed line) and the GRSV2 (dotted line), 
(b) Same as (a) but for the G-S a, G-S b and G-S c parametrizations. 

[2] (a) The unpolarized cross section da/dpr for — 2 < y 1 < 2 as a function for pj, 
at ^/S = 200 GeV as predicted by the GRV (solid line) and CTEQ3M (dashed 
line) parametrizations of the proton distribution functions, (b) Same as (a) but 
for the polarized cross section as given by the GRSV1 (solid line) and GRSV2 
(dashed line) parametrizations of the polarized parton distributions, (c) Same 
as (b) but for the G-S a (solid line) G-S b (dashed line) and G-S c (dot dashed) 
line. The GRSV1 (dotted line) is included for comparison. 

[3] (a) The longitudinal asymmetries for the cross sections displayed in Fig. 2b. (b) 
Same as (a) but for the G-S polarized distributions. The unpolarized distribu- 
tions used is CTEQ3M. 

[4] (a) Comparison of the unpolarized cross section da/dp? as a function of ]P T for 
-2 < y 1 < 2 using the CTEQ3M proton distributions, at \/~S = 200 (dashed line) 
and v^S* = 500 GeV (solid line), (b) The polarized cross section at = 500 
GeV using the G-S a (solid line) and GRSV1 (dot dashed line). Included for 
comparison are the contributions from the qg initiated process. In the case of 
G-S a it is the dashed line and for GRSV1 the dotted line, (c) The asymme- 
try predicted by the G-S a parametrization for the cross section given in (b). 
Included for comparison are the asymmetries as given by the qg (dashed line) 
and qq (dotted line) initiated processes, (d) Same as (c) but for the GRSV1 
distributions. The LO prediction for the asymmetry (dot dashed line) using the 
LO version of the GRSV1 distribution is included. 
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[5] The renormalization/factorization scale (/x) dependence for the unpolarized cross 
section given in Fig.4a, for three different choices of scale, \x = n((p^ 1 ) 2 + 
(pJ 2 ) 2 )/2: 0.5, 1.0, and 2. 

[6] The K-factor for the unpolarized cross section given in Fig. 4a at VS = 500 GeV, 
where K=da(LO)/dp 1 T /da(NLO)/dj? T . 
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